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Air shower simulation using GEANT4 and commodity parallel computing 

L. A. Anchordoqui a *, G. Cooperman b t, V. Grinberg fc t, T. P. McCauley a§ , T. Paul" 1 , S. Reucroft a H, J. 
D. Swain a **, and G. Alversorr 3 ^ 

a Department of Physics, Northeastern University, Boston, MA 02115, USA 
b College of Computer Science, Northeastern University, Boston, MA 02115, USA 

We present an evaluation of a simulated cosmic ray shower, based on GEANT4 and TOP-C, which tracks all the 
particles in the shower. TOP-C (Task Oriented Parallel C) provides a framework for parallel algorithm development 
which makes tractable the problem of following each particle. This method is compared with a simulation program 
which employs the Hillas thinning algorithm. 



1. Introduction 

The steeply falling end of the cosmic ray spec- 
trum now extends up to « 3 x 10 20 eV (see Ref. 
jlj for a recent survey), three orders of magni- 
tude higher than the highest energy achieved by 
hadron colliders. Direct measurements using so- 
phisticated equipment on satellites or high alti- 
tude balloons are limited in detector area and 
in exposure time. Ground-based detectors with 
large apertures make such a low flux detectable 
after a magnification effect in the upper atmo- 
sphere. Namely, the incident cosmic radiation in- 
teracts with atomic nuclei of the air molecules 
and produces extensive air showers which spread 
out over large areas. This indirect method of de- 
tection bears a number of serious difficulties in 
the interpretation of the recorded data. In par- 
ticular, since many variables are involved the pro- 
cesses describing the shower development are in- 
trinsically complicated, numerical simulation of 
the giant cascades has to be performed. The most 
important source of fluctuations in Monte Carlo 
simulations such as CORSIKA || and aires |3| are 
the depth and characteristic of the first few inter- 
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actions, necessarily related to the quality of our 
understandying of hadronic collisions. 

The parts of the shower model which are based 
on electromagnetic or weak interactions can be 
calculated with "good accuracy". The hadronic 
interaction, however, is still subject to large un- 
certainties. In recent years, many models of 
hadronic interactions have been built around ex- 
perimental results, predominantly of pp colliders. 
Extrapolations to higher energies, to small angle 
processes, and to nucleus- nucleus collisions have 
been performed with varying levels of sophistica- 
tion. In particular, the algorithms of SIBYLL (J] 
and QGSJET H are tailored for efficient operation 
to the highest cosmic ray energies.^] The different 
approaches used in these codes to model the un- 
derlying physics can lead to different results when 
applied to the same data 0. Hence, considerable 
systematic uncertainties remain. In addition, a 
significant deviation was recently reported in the 
predictions of aires and CORSIKA, even when 
both programs invoked the same hadronic inter- 
action model to account for the first generation 
of particles M. These discrepancies could result 
from different energy cuts for the hadronic inter- 
action routines, or there may be some discrepan- 
cies in the sampling techniques implemented in 
these programs to reduce the number of gener- 



9 It is important to stress that the latest version of dpmjet 
was improved for operation up to primary energies of 10 21 
eV (per nucleon in he lab. frame) |J . However, it was not 
yet efficiently implemented in AIRES, nor in CORSIKA. 
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Figure 1. Density of air as a function of the ver- 
tical altitude above sea level. 



ated particles. We, therefore, considered it worth- 
while to compare the simulation results of aires 
(which utilizes the Hillas thinning algorithm ||]) 
to that obtained with GEANT4+TOP-C |hJ (top- 
C makes tractable the problem of following each 
particle during the atmospheric cascade). [J] 

2. Air shower simulation 

The most direct way to test the predictions of 
aires and GEANT4+TOP-C is to study the charac- 
teristics of the secondaries generated under sim- 
ilar conditions. Therefore, in order to reduce 
possible differences caused by the implementa- 
tion of hadronic collisions, we analyze the show- 
ers induced by 10 TeV protons. Furthermore, we 
have fixed the energy cuts of GEANT4 according 
to the default values of aires. All shower parti- 

10 GEANT4 is the latest stage in the development of the 
GEANT software, superseding the earlier FORTRAN ver- 
sions with a new object-oriented approach in CH — h 
TOP-C (Task Oriented Parallel C) provides a framework 
for parallel algorithm development [h2|. 



cles with energies above the following thresholds 
were tracked: 750 keV for gammas, 900 keV for 
electrons and positrons, 10 Mev for muons, 60 
MeV for pions and kaons and 120 MeV for nu- 
cleons. The particles were injected at the top 
of the atmosphere (100 km.a.s.l) and the ground 
was located at sea level. In the simulations with 
GEANT4+TOP-C, the atmosphere was defined by 
a stack of 230 layers of increasing thickness and 
decreasing density (with the height above sea 
level). The geomagnetic field was not taken into 
account. The variable density was modeled using 
Linsley's parametrization of the U.S. Standard 
Atmosphere leading to the density profile 

shown in Fig. 1. Each point in the figure cor- 
responds to a layer starting with thickness 50 m 
at sea level and increasing to 1 km at the top of 
the atmosphere. The above description is consis- 
tent with that used in aires. QGSJET was used 
to model the hadronic processes in the simulation 
with aires (we note that at 10 TeV the results ob- 
tained with AIRES+SIBYLL are almost the same). 

In Fig. 2 we show preliminary results of the 
longitudinal development of the total number of 
secondaries that a 10 TeV proton may produce 
after cascading in the atmosphere. We also show 
the longitudinal development of different groups 
of secondary particles (we have considered sepa- 
rately 7, e ± , and ^ ± ). 

At 600 g/cm 2 the ratio between the total num- 
ber of particles produced by GEANT4 and the to- 
tal number of particles produced by AIRES is ~ 2. 
This ratio increases as the shower front gets closer 
to the ground, and roughly reaches the value of 
4 at 800 g/cm 2 . The main reason for this di- 
vergence comes from a difference in the number 
of charged pions produced during the shower de- 
velopment. In the first generation of particles 
(governed by the hadronic interaction) GEANT4 
produces more charged pions than aires. Conse- 
quently, in the first steps, the muon shower profile 
developes faster in GEANT4 than in aires. 

The photopion production threshold is at a 
photon energy of 145 MeV in the proton rest 
frame, with a strongly increasing cross section 
in the region of the A resonances, which decay 
into neutral and charged pions. With increasing 
energy photon-proton collisions are no longer me- 
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diated by resonances, but by the diffractive pro- 
duction of the vector mesons, p and ui, or direct 
multi-pion production. We are currently evalu- 
ating whether the physical models available in 
GEANT4 are sufficient to describe photon-nuclear 
interactions in the energy range of interest. We 
note that GEANT4 provides a framework for intro- 
ducing new physical models if need be. Though 
our (geant4) simulation produces j's above the 
photopion threshold, we observe no charged pions 
from this process. Therefore, as the shower front 
approaches the ground, the neutral pion decay 
channel dominates the cascade, yielding a huge 
number of electrons and positrons. As a result, 
the total number of particles is increased, and the 
position where the shower develops the maximum 
number of particles is shifted to greater depths. 

3. Final Comment 

The above comparison between aires and 
GEANT4 is not yet complete. A detailed discus- 
sion including the geomagnetic field effects, study 
of the 7p-interaction, and the analysis of different 
kinds of primaries at "Auger energies" @ will be 
given in a future publication. 
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Figure 2. Longitudinal profiles of different 
secondary species predicted by aires (o) and 
GEANT4 (★), for showers induced by protons of 
10 TeV. 



